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Abstract —We present a derivation of the radiometer equation 
based on the original references and fundamental statistical 
concepts. We then perform numerical simulations of white 
noise to illustrate the radiometer equation in action. Finally, 
we generate 1/f and 1 / f 2 noise, demonstrate that it is non¬ 
stationary, and use it to simulate the effect of gain fluctuations 
on radiometer performance. 

I. Introduction 

While attempting to understand noise in tube amplifiers, 
Johnson m discovered that “Statistical fluctuation of electric 
charge exists in all conductors, producing random variation 
of potential between the end of the conductor.” The Na¬ 
tional Institute of Standards and Technology (NIST) defines 
a radiometer as a very sensitive receiver, typically with an 
antenna input, that is used to measure radiated electromagnetic 
power 0. Radiometer theory connects Johnson’s and NIST’s 
concepts and leads to the radiometer equation, which is of 
universal importance to radio and terahertz astronomy as 
well as atmospheric studies because it describes the ultimate 
sensitivity of radiometers 0. Unfortunately, most astronomy 
textbook and lecture notes on the radiometer equation do not 
present the details of this relationship, and often provide only 
a cursory overview, which can be misleading depending on the 
background of the reader. In this work, we present a derivation 
of the radiometer equation and explore the effects of non¬ 
stationary noise via numerical simulations. 

II. Derivation of the Radiometer Equation 

A. Variables and constants 

In the definitions below, the typical physical units are given 
in parentheses. 

• A 2 = instantaneous power (watt) = the observed power 
at one statistically independent observation. 

• A = amplitude (Vwatt) = V/VR where V = Voltage 
(volt), and R = resistance (ohm). 

• /3 = equivalent noise bandwidth (Hz) of the radiometer. 
For Gaussian white noise, the equivalent noise bandwidth 
has a uniform power spectral density and represents a 
finite spectral interval of Johnson noise. Bandwidth in 
this description will refer to baseband bandwidth. 

• t — time interval (sec) during which the radiometer 
measures mean power 

• P = mean power (watt). 

• N = number of statistically independent observations. 

• k = Boltzmann’s constant = 1.38*10 -23 (joule/kelvin) 


• T = blackbody temperature of source (kelvin). 

• fi = mean of a distribution. 

• a = standard deviation of a distribution. 

• cr 2 = variance = £(X — p) 2 /N. 

• cr se = standard error of the mean. 

• crp = standard error of mean power. 

• cr T = standard error of blackbody temperature. 

• z = standard normal random variable. A normal random 
variable with p= 0, and a- 1. 

• z 2 = chi-square random variable, one degree of freedom 
which has cr = y/2, and p= 1. 

B. Useful Rule 

If c is a constant, and z 2 a chi-square random variable 
with one degree of freedom, then cz 2 is a gamma distributed 
random variable and the standard deviation is y/2c. A normal 
random variable may be expressed as az + p. 

C. Definition 

A sample is a set of observations of a random variable. Most 
of the statistics involved in this derivation are sample statistics. 
If population statistics are known they may be substituted 
were permissible. The standard error is the standard deviation 
of a sample statistic. The standard error of the mean is: 
cr SE = cr/ VN where cr is the standard deviation of the random 
variable being averaged, and N is the number of statistically 
independent samples used to generate the mean. The standard 
error of the mean is the standard deviation of many sample 
means. 

D. Assumptions 

A is a normal random variable with standard deviation a 
and mean p= 0. The point of reference of this derivation is the 
output terminals of the radiometer antenna or output terminals 
of a resistor. The amplitude A is measured at the output 
terminals. T, A, and P are statistically stationary, meaning that 
the statistical moments and autocorrelation are not functions 
of time. 

E. Derivation 

The radiometer equation is a function that determines the 
standard error of the mean of instantaneous power or temper¬ 
ature. 

• A = az, a statistically independent observation of a nor¬ 
mal random variable with zero mean. Note that A 2 =cr 2 z 2 . 
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. P = ZA 2 /N = cr 2 d. 

• mean power - a 2 — T,(A — p) 2 /N with /i = 0. 

• But also: P = T>a 2 z 2 /N given A = crz + p with p = 0. 

. V2a 2 is the standard deviation of A 2 , given z 2 with p = 

1 and a = y/2. 

• Given that \[2o 2 is the standard deviation of a 2 z 2 the 
standard error of the mean of P is ap = s/2cj 2 /s/N. 
This is the standard error of variance of a normal random 
variable |5). 

• The number of statistically independent samples in a 
signal of bandwidth /3 and in r time is: N = 2 /3r (6). 

• The standard error of the mean power can then be written 
as: 

- cr P = a/2ct 2 /a/]V 

- CTp = \[2cr 2 / a/2/3t 

- crp = G 2 !\fWr 

• Given Johnson’s m formula P = kT/3, the standard error 
of temperature is: <tt = T/yf/Pr. 

III. Numerical simulation of the radiometer 

EQUATION 

We have used the numerical python packages (numpy and 
scipy) to simulate and analyze a radio frequency signal as 
a zero mean white noise datastream of 10 million points. 
We set the standard deviation of A = 2 in arbitrary units, 
resulting in a power of 4. In Figure [I] the left column shows 
the amplitude (voltage) signal vs. time, a histogram of its 
probability distribution function (PDF), and a histogram of 
its normalized PDF. The middle column shows the same 
quantities for the amplitude squared (i.e. power) signal. In 
this case, the power PDF is described by a Gamma function 
(shown by the red line), while the normalized power PDF is 
a chi-squared function with 1 degree of freedom. In the right¬ 
most column, we show PDFs of the observed power which 
simulate a series of observations of increasing length. In the 
top right panel, we break the datastream into 10 5 samples 
each containing 100 observations. With only N=100 samples 
in each observation, the width of the PDF is large, showing 
that the signal to oise ratio (SNR) of one observation is only 
7.07. Increasing the number of samples to 7V=1000 and then 
10000 shows the increasing precision on the measurement of 
the power level to a SNR of 22.3 then 70.2, consistent with 
the prediction of the radiometer equation. 

IV. Numerical simulations with non-stationary 
noise 

A. Importance of stationary vs. non-stationary noise 

The radiometer equation is based, in part, on the assumption 
of statistically stationary mean and variance of the Gaussian 
signal from the thermal source of radiation observed by the 
radiometer. This provides a good estimate of mean power and 
temperature. The estimate fails, at some level of precision, 
due to the time varying gain and noise temperature of the 
radiometer. Simple solutions to this issue can be found if 
the time variation is predicable, for instance a linear trend 
like ones age, or a cyclostationary variation, like the average 
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Fig. 1. Simulations of the radiometer equation with a stationary white noise 
signal, as described in section §[In] The left column shows the amplitude time 
series along with its mean, standard deviation and PDFs. The middle column 
shows the same things for the corresponding power signal. The right column 
shows how the PDF of the measured power becomes narrower around the true 
power (marked by the vertical green line) as the integration time increases, 
consistent with the prediction of the radiometer equation. 

temperature in an hour of a day, which shows both diurnal and 
annual periodicity. Unfortunately the variation of radiometer 
gain and noise temperature are random. 

A common model for this behavior is a random variable 
whose power spectral density has a power law dependence. 
Johnson had to take a baseline measurement of the amplifier 
response with no signal at the input, and subtract this “zero 
deflection” from his measures of thermal noise. Dicke devel¬ 
oped the Dicke switch method of radiometer gain calibration 
to overcome the gain variation m . Power law noise is typically 
found at the low frequency end of the power spectral density 
of an amplifier’s output power m The power is a function of 
frequency and may be represented by P oc f~ a where a is 
slope of the power spectral density. 

B. Generating 1/f noise 

To simulate 1 / f a noise with different values of a , we use 
the C code of Paul Bourke Q which creates noise data by 
the “fractional Brownian motion (fBm)” method. The fBm 
method involves creating frequency components which have a 
magnitude that is generated from a Gaussian white process and 
scaled by the appropriate power of a (the phase is uniformly 
distributed on [0, 2tt]) and then inverse Fourier transformed 
back into the time domain. Figure 3 shows that the power 
spectrum for the resulting a = 1 is indeed 1/f. 

One feature of noise with a > 0 is that the variance 
diverges with time, in contrast to white noise whose variance 
approaches a constant value (Figure 4) 

In subsequent figures, we show how this noise differs from 
Gaussian white noise in a number of different representations. 
In Figure 5, we show the structure of a portion the time-series, 
which exhibits a quasi-periodicity. 

1 See http://paulbourke.net/fractals/noise 
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Fig. 2. The left column shows a stationary datastream the mean and variance 
do not change with time. The right column shows a non-stationary datastream 
in which 100 consecutive samples can sometimes be all positive (first row), 
or all negative (second row) or be centered at zero (third row). Although both 
datastreams produce a Gaussian PDF when viewed in aggregate (bottom row), 
the variance diverges with time when a > 0, as shown in Fig. [4] 
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Fig. 5. Plots of ten thousand power law values vs. time for a signal with 
a = 1 (left panel) and a = 2 (right panel) show a periodicity which arises 
from a low frequency component of the noise. If we were to zoom into the 
plot, an image with similar features, though decreasing amplitudes would 
appear. 
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Fig. 3. Power spectrum of white noise (a = 0, left panel) compared to a = 1 
noise (right panel) The red line has a slope of zero and -1, respectively. 



Fig. 4. Variance of a signal vs. the logarithm of the total time examined 
for three different random signals: a = 0 (solid line), a = 1 (dashed line), 
a = 2 (dotted line). Signals with a > 0 diverge with time. 


The auto-correlation functions (ACFs) are compared in 
Figure 6. As expected, the ACF of white noise has signal 
only in the zero lag, and is zero at all other frequencies. In 
costrast, a power law random variable displays a decreasing 


Fig. 6. Graphs of the autocorrelation function (ACF) of three random 
variables. Left panel: white noise (a = 0); Center panel: power law noise 
with a = 1; right panel: power law noise with a = 2. Different time samples 
are shown in different colors. Note that the ACF changes with time for a > 0. 

correlation with lag, and the details change with time. In other 
words, the ACF changes when different independent samples 
are examined. 

In Figure 7, we show the Allan variance of a signal which is 
primarily white noise but has a portion of a > 0 noise added 
to it. In Figure 8, we show how the Quantile-Quantile (QQ) 
relation for the a = 1 and a = 2 signals compares to that of 
a uniform distribution and a Cauchy distribution. 

Gain fluctuations are an unfortunate but common feature 
of low noise amplifiers such as HEMTs 0, ED- Here we 
simulate the effect small gain fluctuations by imposing a small 
fraction of power law noise with a = 1 onto a white noise 
radiometer signal. The graph of frequency distributions of the 
power law random variable over time shows that estimates 
of the mean value have a greater dispersion than would be 
expected from the standard error of the mean. The effect 
increases as a greater amount of gain fluctuation is inserted. 
The S/N ratio no longer improves as fast as 1/Vtime. 

V. Further reading 

A discussion of the radiometer equation in terms of single 
dish radio telescopes is given in the summer school lectures 
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Fig. 7. Allan variance of a signal composed of the sum of white noise plus 
some a = 2 noise. 
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Fig. 8. The probability density function of the power law noise in the 
models that we use appear similar to Gaussian (normal) for values of below 
2 and more like a uniform probability density function for above 2. This 
behavior can be displayed in the graph of the normal Quantile-Quantile (QQ) 
plot (upper panels), where a normal distribution would lie perfectly on the 
diagonal. This can be compared to the QQ plots of the uniform distribution 
(lower left) and Cauchy distribution (lower right). 


GD Similarly, the application of the radiometer equation to 
cross-correlating interferometers is treated by several authors 

G2, 113, £E3. 


VI. Future work 

Our future work on this topic includes simulating the role 
of the radiometer equation in a two element interferometer, 
and examining the statistics of the digital data recorded by 
the Submillimeter Array (SMA) fl5l . lH6l in Hawaii. 
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